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Based on the first principle calculation, a Lagrangian for the system describing quarks, 
gluons, and their interactions, is constructed. Ascribed to the existence of dissipative 
behavior as a consequence of strong interaction within quark-gluon plasma (QGP) mat¬ 
ter, auxiliary terms describing viscosities are constituted into the Lagrangian. Through 
a ’’kind” of phase transition, gluon field is redefined as a scalar field with four-vector 
velocity inherently attached. Then, the Lagrangian is elaborated further to produce the 
energy-momentum tensor of dissipative fluid-like system and the equation of motion 
(EOM). By imposing the law of energy and momentum conservation, the values of shear 
and bulk viscosities are analytically calculated. Our result shows that, at the energy level 
close to hadronization, the bulk viscosity is bigger than shear viscosity. By making use 
of the conjectured values g/s l/47r and C,/s ~ 1, the ratio of bulk to shear viscosity is 
found to be I^/t] > 47r. 
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1. Introduction 

Relativistic heavy ion collision experiments have confirmed that a strongly coupled 
quark-gluon plasma (QGP) has been reproduced. Lots of new phenomena from this 
new state of matter, like its small transport coefficients and the real phase transition 
mechanism from deconfined quarks to confined hadron, have currently become hot 
research topics. 
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Strong interaction typically indicates the existence of small transport coefficients 
of a system in order to ensure the efficiency of momentum transmission. In such kind 
of system, shear and bulk viscosities are naturally become its inherent properties. 
Understanding these properties firmly has been considered as a prerequisite knowl¬ 
edge for dynamical description of plasmas, astrophysics, early universe cosmology, 
and other branches of physics. 

The dynamics of quark-gluon plasma that intuitively been considered to live in 
the realm of quantum chromodynamics (QCD) has been studied through the lattice 
theory.®^ Nevertheless, some simplifications are always required in order to reduce 
the calculation complexities as well as to cope with the limited computing capacity 
of the current state-of-the-art computer. On the other hand, experiments of relativis¬ 
tic heavy-ion collision also exhibit a collective flow behavior. The typical character 
which belongs to fluid dynamics is considered as one of the validating aspect for 
QGP to be modeled as a relativistic hydrodynamics system.l^lti^In the development, 
the discovery of new phenomena from experiment such as the shock wave in form 
of Mach cone has guided researchers to materialize the gluonic-plasma model.^^i^ti^ 
Recent achievements in hydrodynamics description of relativistic heavy-ion colli¬ 
sion have inspired more researchers to use relativistic fluid dynamics approach on 
searching for viscosities of QGP. 

There is also a hybrid model that describes the relativistic hot plasma by uni¬ 
fying electromagnetic fields with flow field. The unification is represented by the 
effective field strength tensor, -|- {rn/q)S^v, combining appropriately the 

electromagnetic and fluid fields.l^ 

In this paper, we follow this hybrid concept, except that all fields are constructed 
by following the first principle. In this framework, the strongly interacting system 
of QGD is considered as a macroscopic fluid system rather than the result of many- 
particle interaction. This scenario is, somehow, inspired by the works of Jackiw 
et where they took an analogy of a non-commutative quantum system in 

developing a mathematical description of fluid structure, that is still consistent 
when applied to the microscopic fluid system (non-Abelian fluid). Instead of taking 
the same route, we start with different way, i.e., by following the first principle. 
To this end, a viscous fluid-like QCD Lagrangian is constructed. By redefining the 
gluon fields, a sort of phase transformation on the Lagrangian appears. We then 
obtain an expression that can be interpreted as a microscopic interacting fluidic 
system. This strategy also produces a hybridization that combines the flow field 
and the charge field in one system. Further, by making use of the Lagrangian, the 
formulation for bulk and shear viscosities can be obtained. At the energy level near 
to the phase transition, the result of our calculation exhibits a large ratio of bulk 
to shear viscosity, just as what is being seen at RHIC. This result corroborates the 
new findings on QCD plasma. 

Our paper is organized as follows. In Sect. [2] we briefly review the basic theory 
of fluid QCD, and explain the phase transition from the deconflned particle field to 
the flow field. In Sect. (H the shear and bulk viscosities of gluon dominated QGP 
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are formulated as a function of the energy density. The obtained results will be 
discussed in Sect. |4l Finally, we will summarize our finding in Sect. [5] Throughout 
this work we use the natural units, i.e., = c = 1. 

2. Theoretical background: The gluon dominated QGP 

In general, the Lagrangian of the viscous QGP can be assumed to consist of a 
nondissipative (ideal) and a dissipative part. Mathematically, we can write 

^ — .^ideal “t“ .^diss. t (1) 

where the ideal part reads 

Adeai = iQl^d^Q -mgQQ- + 9gJg^A^^- (2) 

Here Q and Q represent the quark and anti-quark triplet, 7 ^ is Dirac matrices, and 
mq is the quark mass. The kinetic term for gauge boson is built inside its field 
strength tensor, i.e., —The same role is taken for the gauge boson 

I/“ of SU(3) gauge group, the kinetic terms — di,U^ + 

where gp is the strong coupling constant. The last two terms, = QTpjfj,Q and 
Jgii = Q7fj,Q, represent the quark currents from the SU(3) and U{1) gauge groups, 
respectively, and go denotes the coupling constant from U{1) gauge group. 

Being a Lagrangian of QGD, Tideai is constructed from the fact that the derived 
energy-momentum tensor has an identical form as the energy-momentum tensor of 
an ideal fluid, as shown by Ref. 18. It is the dominant part of the Lagrangian that 
describes the unification of fermions and bosons from different gauge groups, that 
preserves the SU(3 )f G U(1)g gauge symmetry. It represents a non-Abelian SU(3) 
fluid fields containing quarks and anti-quarks that interact with an electromagnetic 
field. 

The nonideal or dissipative part of Eq. o can be written as, 

Tdiss. = -h - U‘^''U°-^deUa^) 

-T^iC - . (3) 

This part is constructed from the relativistic hydrodynamics viscous energy- 
momentum tensor developed by Landau and others.^^In Eq. T] and C denote 
the shear and bulk viscosities, respectively, whereas T“ is the structure constant of 
the SU(3) gauge group. 

Therefore, the complete Lagrangian can be written as 

£ = iQrd^Q - mqQQ - 

-b d''u^ - - u^'^u^^^deUa^) 

-T^iC - lv)deU^%l - . (4) 

By introducing the viscous terms into the Lagrangian, it is corollary representing 
a ’’near to ideal” fluid, while the system itself is near to the equilibrium. Gauge 
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symmetry of the Lagrangian is preserved until the dissipative terms are added. 
This is understandable as the fact that gauge invariance is broken when the fluid¬ 
like system leaves its ideal condition. 

To introduce the flow property into the fluid field, a construction of relativistic 
form is in order. It has been argued that [/“ should be rewritten in a particular 
form of the relativistic velocity as C/“ = (C/q,U“) = Here, = 7 ( 1 , v) 

and 7 = (1 — with tjf' = (jf'ix) represents a dimension one scalar field. 

Equation of motion (EOM) for a single gluon field from the above Lagrangian leads 
to a general relativistic fluid equation.^l^ This fact explains that a single gluonic 
field t/“ may behave as a fluid at certain phase and as a point particle at hadronic 
state with a polarization vector in the conventional form of C/“ = This is 

considered as a kind of ’’phase transition” between fluidic and hadronic states. 

hadronic state QGP state . 

'---" S-,-^ 

A more detailed discussion on this topic can be found in Ref. 21. 

By applying the Euler-Lagrange equation to Eq. o, and ignoring the minor 
contribution from the dissipation part, we obtain 

^(7v</>“) + V(7r) = -9f j dx(J“o + F“). (5) 

In the non-relativistic limit, i.e., 7 ~ l-|-^|vp and </>“ 1, for a constant Eq. dSD 

reduces to the classical equation of motion of fluid dynamics, i.e., 

9v f 

— -b (v • V)v = -gp j) dyi{JpQ + Fq) |„on-rel, (6) 

where JpQ is a covariant current induced by the existing quarks surrounded by and 
interacting with gluon ’’fluid”, while is induced by the fluid self-interaction and 
the interacting gauge fields A“. This shows that from a certain point of view the 
Lagrangian describes a general relativistic fluid system interacting with another 
gauge field and matter inside. Besides that it also fulfills the basic presupposed 
idea that in comparison with the size of a system, the relatively long-distance, 
low frequency behavior of any interacting theory at finite temperature could be 
described by the theory of fluid mechanics 

From Eq. m it is obvious that the flow characteristic describing a plasma is 
represented by the gluon fields and, accordingly, the system that we are working 
with is a gluon dominated QGP. Therefore, all nongluon terms or terms that have 
no interaction with gluon can be omitted. As a consequence, £diss does not change 
and £ideai is simplified. Now the Lagrangian reads 

c = 

+ d''U^- - U‘^''U°-^deUa^) 


( 7 ) 
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Here £ideai = Besides describing the kinetics of gluons, 

it also indicates the self-interaction between small number of quark and anti-quark, 
the self-interaction between gluons, and the interaction between quark with gluon 
” fluid”, where it is immersed in. It should be noted that the quarks and anti-quarks 
feel the electromagnetic force due to the U{1) field but the size is suppressed 
by a factor of ejg = \Jajus ^ 0(10”^). 


3. Shear and bulk viscosities 

From the Lagrangian one can also derive the energy momentum tensor density by 
means of Tfiv = ■ For the ideal term £ideai, one immediately geti^ 

^deai = + 2gprp^U; 

= [2gpg^,rp^U-^^ + glr^^r'^^UlU^U'^PUt] 

Jp^U'^^ in Eq.® can be obtained by considering the EOM (Dirac equation) of a 
single colored quark [Q) or anti-quark (Q) with 4-momentum . Since the solution 
of the EOM is Q{p,x) = q(jp) exp(—ip • x), one immediately gets q'y^q = 4p^. 
Assuming that all colored quarks/anti-quarks have the same momenta and the 
velocity of gluons and quarks inside the gluon sea are homogeneous, approximately 
Jp^U°'^ oc Ap^T°'U^ = AmQT°-(f>°- since Ufj,u^ = = 1. Then the ideal part of 

energy momentum tensor in the function of field reads, 

^deal = [SgpruQT-r + 

-[AgpruQT^r - • (9) 

After obtaining 7)ji/ ideal that represents a perfect fluid, we can derive the dissipative 
part of energy-momentum tensor 7)ii/ diss from Ldiss, i-e., 

d^ss = -vT‘^id,u^ + d^u: - - u^u^^deU:) 

+|ryr“5eC/“^(p^, - c/;c/“) - - c/;c/“). ( 10 ) 

The structure of the outcomes reveals a typical form of the relativistic hydrody¬ 
namics viscous energy-momentum tensor.^l^ Further, we assume that the gluon 
color states are homogeneous, i.e., (j)°- = (f> for all a = 1, • • • , 8 . As the follow-on to 
this assumption, the generator r“ can be compactly written as T = 
also the 2nd and 4th terms of Tpi, ideal can be omitted due to the completely anti- 
commutative property of the structure constant Now, the complete energy 
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momentum tensor that consist of ideal and dissipative part reads, 



■qT{dy{u^4i) +d^{Ut,(j)) - (j))de{u^4')^ 


+ ‘^'nTde(u^4>)gfj,^{l - (f^) - C,Td^{u^(t))g^^{l - 0^) 


( 11 ) 


The total energy momentum tensor is obtained by integrating Eq. m over the 
total volume of the space-time under consideration. It reveals the collective gluons 
flow in the system. 

The shear and bulk viscosities can be obtained by imposing the energy and 
momentum conservation to the energy-momentum tensor. The conservation of en¬ 
ergy is expressed as 9oT°° -I- dkT^^ = 0.1^ After some algebra and assuming that 
= 0 at a very short traveling distance and very short existence time of QGP, 
the conservation of energy can be formulated as 



( 12 ) 


where 


c° = SgpmQ {do(j))u°u° - i(9o</')ff°° + {dk(t>)u°u'") , 


a° = 29 o5°(w°(/)) — 2{yPu^{do4>^){de(j)) -I- n°u°u®0^(9o9e^i)) 

-^(5o9e(u^(/)))g°°(l - (^^) -h ‘^{de{u°‘(j)))g°°{dQ(t)'^) + dkd'^{u°4>) 
+dkd^{u^4>) - u'"u^{dk(j)'^)de{u°(j)) - u'^u^4>^{dkde{u°(t>)) 

-vPu^{dk(l>^){de{u’^(t>)) - («'=(/>)), 


{dode{u^(j)))g°°{l - - {de{u^(j)))g°°{do(j)^). 


The conservation of momentum, -I- -I- dkT^^ = 0, can be re-expressed in 

a simpler form as 



( 13 ) 


where 




z* = SgFrriQ [{do(j))u^u^ + {di(j))u^u^ 
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= 5oa°(wV) +5o5*(w» 

—u°u’^{do(j)^)ds{u^(j)) — vPu^(jP'dQds{u'‘4>) 

— — U^u‘^(j)^dode{u°(j)) 

+2did’'{u''4>) — 2M®M®(9i(/)^)i9e(u*(^) — 2u''u^(j)^didg{u’‘(j>) 

+i9fei9*(u^(/>) — u^u^{dk4>^)de{u^(j)) ~ dkde{u'‘4>) 

-u'u’'{dk'f)de{u’^(l)) - u^u^'ifdkdeiu^^cj)), 

- </>") - 


With the availability of Eqs. (HU and dH, the bulk and shear viscosity can be 
expressed explicitly, i.e.. 


c 


c°a;* - a°z* 


(14) 


a°y^ - x^b° ■ ^ ^ 

The expression of the scalar field (j) can be obtained by solving Eq. ([S]). For 
simplicity, but without losing generality, the partial differential equation is assigned 
to live in two-dimension, (j) = 4>{t,x), and the gauge bosons from each gauge group 
are assumed to be homogeneous. As a consequence, the differential equation becomes 

dt{-fV(j)) + da,{-i4>) = -gFiJpo + Fo), (16) 

with dt = d/dt, dtx = jdtdx, gp and go are the gauge group SU(3) and U(l) 
coupling constants, respectively. Here, JpQ = Q^^T'^Q = and 

Fq = (jidtfj) - {gG/gF)A-fdt(j) + - {go/gF)A'^dx(l>. (17) 

Thus, Eq. (IT6l) can be written as 

bdtx(l) + ldxx(l> + C + I34>dt4' + I34'dx4> - - \dx4> = 0, (18) 

with S = 7|v|, ^ = gFQ"/fiT°‘Q , P = *7^|vp, and A = {gc/gF)A^. This equation 
can be solved analytically when it is converted to an ordinary differential equation 
form by redefining </> = (j){x — iEt), where x — iEt = z. Then, it follows that dx = dz, 
dt = —iEdz, dtx = —iEdzz, and dxx = dzz- Therefore, Eq. (IT^ becomes, 

Otlpzz - a2(j)4'z + Oi^pz = i, (19) 

with ai = —iE6 -|- 7, a2 = iEP — P, and = iEX — A 

Since the particular solution contributes only to transient condition, here we try 
to obtain an analytic homogeneous solution. 


aipzz - ot2ppz + OL^Pz = 0, 


(20) 
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Fig. 1. Bulk viscosity of the gluon dominated QGP as a function of the temperature calculated 

from Eq. Gl¬ 


and by defining 


(j) = 



we obtain 


aiUJzzz - \a2{ujl)z + OfiUJzz = 0 . 
Equation m can be integrated to obtain 


dujz 


J {\a20jl - asujz) 

and by making use of the standard integral formula 

f dx ln[(6a; — a)/a;] 


dz 


6x2 


we arrive at 


• = UJz = 




Q;2eGia3 _ 2 ^C)Lzz/Cii'\-C'2.C)Lz 


( 21 ) 


( 22 ) 


(23) 


Eor the sake of simplicity, the integration constants Ci and C 2 are set to zero, which 
reduces the field to 

, 2a2 


0-2 — 2e“3(^“*-®0/ai ■ 


(24) 
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Fig. 2. Shear ( 77 ) viscosity of the gluon dominated QGP as a function of the temperature calcu¬ 
lated from Eq. GUI. 


4. Results and Discussion 

To visualize our result for the gluon dominated QGP at hadronization temperature, 
we take the temperature range between 150 and 250 MeV. Note that 150 MeV is 
the hadronization temperature that commonly adopted in the relativistic hydrody¬ 
namics description of QGP, while the values above 200 MeV are usually considered 
in the lattice QCD computation for QCD plasma. The other variable values are also 
chosen with respect to the condition where gluon dominated plasma is conjectured 
to exist. The particle velocity v is assumed to be near the velocity of light, whereas 
At = 1 X 10“^^ s. The space variable x is taken from 10“^ to 10“^ nm and the 
quark mass mg = 0.5 GeV. The energy of gauge field A is set to 100 GeV. This 
value is taken with an assumption that the wave length of gauge field is about two 
order smaller than the radius of quark rg, where rg ~ 0(10“^®) m. Then, by using 
go ~ 1/137 and gp ~ 1.00, the ratio between weak coupling constant and strong 
coupling constant gc/gp is found to be 0.0072. It should be noted that from all of 
the variable values where gluonic plasma is conjectured to exist, the values of quark 
mass mg and the energy of gauge bosons A span a large range. However the result 
in the figures does not much depend on the values of these variables. The reason 
is that the value of mg and A can be reciprocally altered due to the assumption 
that the higher the energy level of gauge boson in a circumstance, the smaller the 
possibility of heavy particles to exist. If A is set to a bigger value, then mg has to 
be set to a smaller one, and vice versa. 

Figure [T] shows that the obtained value for bulk viscosity increase, especially at 
T very close to T ~ 0.15 GeV. It is not totally similar with the result which is 
obtained from the calculation that utilizing quantum kinetic theory, where near to 
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Tc the bulk viscosity rises asymptotically. 

Figure exhibits the calculated shear viscosity, which has a small increment 
near the hadronization energy level. Compared to the bulk viscosity depicted in 
Fig.m the shear viscosity turns out to be one order smaller than the bulk viscosity. 
Qualitatively one can see that ^ ~ C>{10), and r] ^ 0(10°). With the conjectured 
value of shear viscosity to entropy ratio , rj/s > (Kovtun-Son-Starinets 

bound in natural units) one can conclude that the entropy of QGP is also of the 
order of s ~ 0(10), and C/s > 0(10°). This result is comparable to that of Tuchin 
et. al. for the QCD plasma, where C/s ipHHH -vvhich indicates that for a hot 
QCD matter there is a significant increase of the bulk viscosity at the vicinity of 
critical temperature, where in this case the adopted critical temperature is about 
0.28 GeV. It might be seen as an indication that at the hadronization process, bulk 
viscosity plays a more important role than shear viscosity. At least, mechanically, 
this large C brings the tendency that the fluidic QGP becomes unstable.^ 

Gomparison between Fig. [1] and Fig. [2] reveals a linear relation between C and 
t] at about hadronization energy level. It can be simply written as C ~ c ry, where 
c ^ 50. Equation ([9]) shows that in the gluon fluidic system there exists no self¬ 
interaction between gluon flleds. Instead, the interaction originates from the matter 
glnon relation. The large value of C originates from the strong coupling constant in 
the gluon dominated system. This is just the opposite to the huge ratio of rj/C, that 
appears at weakly coupled system.l^ This is also related to the fact that for any 
nonconformal field theory the bulk viscosity will always emerge and a non-Abelian 
gauge theory like the gluon dominated QGP that is modeled in this discussion 
becomes nonconformal at hadronization energy level.l^ 


5. Summary 

Based on the first principle, an analysis started from the construction of a viscous 
fluid-like QGD Lagrangian has been performed. The phase transition between stable 
hadronic state and highly energized QGP state has been investigated. It is argued 
that inside the hadronic state the gluons behave as point particles and their prop¬ 
erties are determined by their polarization vectors. However, in the hot QGP state 
the gluons should behave as fluid particles and characterized by their relativistic 
velocities and their inherently attached scalar fields. By investigating the energy- 
momentum tensor for viscous fluid particles, a detailed derivation of the bulk and 
shear viscosities close to the deconfinement temperature has been performed. The 
result of this investigation corroborates the findings obtained in the previous studies 
that make use of the QGD low-energy-theorem approach. 

Another interesting issue related to the result of the present analysis is the 
calculation of cross-section. In our approach presented in this paper, which is based 
on hydrodynamics, calculation of the cross-section seems to be difficult. In principle, 
it would be possible to relate the result obtained for the viscosity with the prediction 
of kinetic theories, where the transport coefficients like bulk and shear viscosities in 
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a gluon dominated QGP are connected to the gluon cross-sectionP^ Nevertheless, 
this would require a much more comprehensive and careful investigation, which is 
planed for the future works. 
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